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Abstract. We examine the prospects of discrete quantum walks (QWs) with trapped ions. In 
particular, we analyze in detail the limitations of the protocol of Travaglione and Milburn (PRA 2002) 
that has been implemented by several experimental groups in recent years. Based on the first realization 
in our group (PRL 2009), we investigate the consequences of leaving the scope of the approximations 
originally made, such as the Lamb-Dicke approximation. We explain the consequential deviations from 
the idealized QW for different experimental realizations and an increasing number of steps by taking 
into account higher-order terms of the quantum evolution. It turns out that these become dominant 
after a few steps already, which is confirmed by experimental results and is currently limiting the 
scalability of this approach. Finally, we propose a new scheme using short laser pulses, derived from 
a protocol from the field of quantum computation. We show that the new scheme is not subject to 
the above-mentioned restrictions, and analytically and numerically evaluate its limitations, based on 
a realistic implementation with our specific setup. Implementing the protocol with state-of-the-art 
techniques should allow for substantially increasing the number of steps to 100 and beyond and should 
be extendable to higher-dimensional QWs. 
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1. Introduction 

Random walks are powerful models that allow to describe, understand and make use of 
stochastic processes occuring in a wide variety of areas [BIS]- Models describing related 
processes in the quantum world are called quantum walks (QWs) [3]. 

An implementation of a discrete (-time and -space) random walk on a line requires 
two basic operations. The coin operation, with the random outcome of Heads or Tails, 
is followed by the shift operation to the left or right, depending on the outcome of 
the coin toss. After steps the walker will therefore have followed randomly one 
out of many possible paths, with the probability for its location being given by a 
binomial distribution centered around the starting point. The average displacement 
of the walker, i.e. the standard deviation of that distribution, increases with the square 
root of A^. The quantum mechanical version replaces the probabilistic coin toss by 
a deterministic operation. It prepares the quantum coin in an (equal) superposition 
of Heads and Tails. As a consequence, the walker performs the conditional step in 
both directions simultaneously. The walker follows all paths during this deterministic 
(and thus reversible) process, allowing for constructive and destructive interferences at 
subsequent crossings. The probability distribution of the position of the walker is due 
to these interferences substantially different from a binomial one. In particular, the 
average displacement of the walker scales linearly in A^. 

Quantum walks have been thoroughly investigated theoretically and several 
applications for QWs have been proposed, for example in terms of quantum computing 
[1]. Many classical algorithms in computer science make use of random walks for 
sampling purposes. Algorithms of that kind might get substantially speeded up by 
quantum versions of the random walk, where all possible paths are tested in parallel, 
potentially providing a similar gain as the prominent example of Grover's search 
algorithm [5]. In addition, QWs can be interpreted as the one-particle sector of a 
quantum cellular automaton, which is a fundamental model of a quantum computer [6J. 
Furthermore it has been shown that QWs themselves are suitable for universal quantum 
computation [7] and different aspects of quantum information processing [H |9l [TO] . 

In a different context, QWs can be exploited as prototype models for intriguing 
transport processes in nature. One examples is the energy transfer in photosynthesis 
with an efficiency of close to 100% pT| IT2] . a performance that is not achievable 
classically. Other examples are the creation of molecules in interacting QWs [I3j and 
effects like Anderson localization and diffusive scaling in disordered QWs [T11[T5]. Here, 
QWs might be suited for experimental quantum simulations to provide deeper insight 
into complex quantum dynamics. Additionally, even relativistic effects can be considered 

Promising attempts at their implementation have been performed for the discrete 
and the continuous versions of QWs. Important aspects of QWs have been realized in 
a nuclear magnetic resonance experiment [17] using the internal degrees of freedom of 
molecules to span the coin and position space. An implementation based on neutral 
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atoms in an optical lattice [THl Hn] has resulted in an experiment [20] where the lattice 
sites in a standing wave of light span the position space of the walker/atom, two 
electronic states encode the two coin states and a state-dependent optical force provides 
the conditional shift. Other proposals consider an array of microtraps illuminated by 
a set of microlenses |21], Bose-Einstein condensates [22] and atoms in cavities |23] . 
Photons have mimicked single walkers on the longitudinal modes of a linear optical 
resonator [23] and in a loop of a split optical fibre [21]. Single [25] and two time- 
correlated photons [261 EZj have recently been travelling and interfering in a lattice of 
optical waveguides. Travaglione and Milburn [28] proposed a scheme for trapped ions to 
transfer the high operational fidelities [29] obtained in quantum information processing 
into the field of QWs. While coin states and steps are operated similar to the atoms in 
the optical lattice, the position is encoded in the motional degree of freedom of the ion(s), 
which oscillate in a quantized harmonic trapping potential. The proof-of-principle has 
been performed by our group by the implementation of a discrete, asymmetric QW 
of one trapped ion along a line in phase space [30]. Recently, the proposal has been 
theoretically refined [31] and experimentally extended to an increased number of steps 

m- 

All of the above-described systems and their related protocols of implementation 
are severely limited in the total number of steps due to a lack of operation fidelities or 
even fundamental restrictions. However, a larger number of precisely performed steps 
is the crucial prerequisite to exploit QWs for the envisioned applications. For the case 
of trapped ions, the limit of coherent displacements to states inside the Lamb-Dicke 
regime has already been foreseen [28], experimentally observed in a different context 
[33l l34] , and confirmed by us [35] and others [32] . 

In this paper we substantially extend the description of the experimental 
implementation of the asymmetric QW with three steps [35]. We carefully analyze 
the effects that arise when approaching the fundamental limitations of the proposed 
protocol [28] after the related, severely restricted number of steps for different step 
sizes. We consider higher-order terms to the soon overstrained approximation building 
on the work of p33]. In parallel, we experimentally confirm the essentials by further 
investigating our results [35j , which already lead into a regime where the refined theory 
is required. Finally, we develop a novel protocol for a QW, based on a scheme from 
the field of quantum information processing using photon kicks [36l [37], to overcome 
these restrictions and to allow in principle for hundred(s) of steps, extendable to QWs 
in higher dimensions. 

The paper is structured as follows: In section |2] we give a theoretical description 
of the QW as it has been realized in our experiment, similar to the original proposal 
[28] , and analyze issues concerning non-orthogonality of the position states. In section 
Owe describe the experimental method of realizing the necessary operations of the QW 
and analyze the limitation of the position space to the LDR of the optical dipole force. 
In section 13.51 we describe the experimental methods of the ion state detection. In 
section H] we describe the experimental procedure, in particular the determination of 
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the relevant parameters. In section E] we summarize the results and the limitations of 
the implementation of the QW. Finally, in section |6l we propose the implementation 
of the shift operator with short laser pulses (photon kicks) and the extension to higher 
dimensions. 



2. Theoretical considerations 



In the following we give a theoretical description of the discrete QW on a line as it is 
realized in our proof-of-principle experiment for the first three steps. 
Consider a Hilbert space vector 

\^) eM = m,orn®^moUon. (l) 

Hcoin denotes the coin space with basis states 



H)- : i m= : , (2) 




encoding the coin states, Heads and Tails. MmoUon is the infinite dimensional phase 
space, related to a harmonic oscillator. We encode the discrete positions as coherent 
states 

CO ^ 

K)=e-l-lV^5:^|n), (3) 



where G Z and at = k ■ Aa with Aa G C The states |?t.) denote the (orthonormal) 
Fock states. For the QW the distance |Aq;| between neighbouring positions in phase 
space is of importance, whereas the argument of the complex number Aa can be chosen 
to be constant for all steps and therefore is irrelevant. 

Concerning the notation of the position states, we will use the following convention. 
Ideally, the position states are coherent motional states as described above. In 

the experiment the position states will contain a small amount of motional squeezing 
(Sect. 14.11) . These states will be denoted as |5fc). Additionally, whenever necessary, 
we will distinguish between position states generated in a numerical simulation, |5f), 
and experimentally, |5f). Further, we will generally use the superindices ^ and ^ to 
distinguish between simulation and experiment, whenever necessary. For the simulation 
we will use 3SB (See sect. 13.4.21) . if not stated differently. 

Ideally, the initial state of the QW is chosen to be 

|7/>o) = |T) ® |ao = 0) . (4) 

Each step of the QW is described by the subsequent application of the coin operator C 
and the shift operator S. Thus the state after steps is given by \iIjn) = {S- C)^ li^o)- 
The coin operator C is defined as 

/TT \ if 1 e^'^\ 
C = i? 0j = -= ^ ® (5) 



Experimental Simulation and Limitations of QWs with Trapped Ions 



5 



position index 

-3-2-10 1 2 3 





-3-2-1 0,1 2 3 

position index 



Figure 1. The asymmetric QW on a line, (a) Schematic of the QW. Starting at the position with 
index 0, the walker tosses a coin and does a shift to the position index 1 (—1) for the coin showing Tails 
(T) (Heads (H)). If the coin is in a superposition of T and H, the walker shifts into both directions 
simultaneously, taking the part related to T to position 1 and the one related to H to position — 1. 
Further coin tosses from T [H] cause superposition states with different relative phases, that is, H + T 
and H — T respectively. The relative phases allow for interferences of the wave function between 
different paths. The first interference occurs during the third step, taking the walker from the state "02 
to tp3. At position the coin toss results in constructive (destructive) interference for T{H), illustrated 
by the bold (dashed) arrow for the subsequent shift, (b) We implement the QW with a trapped ion, 
where the position space is encoded into the co-rotating phase space (Re{a), Im(Q;)) of the axial 
normal mode of motion. The positions k are represented by coherent states \ak), which are aligned 
along a line in the co-rotating phase space. Two electronic (hyperfine) states of the ion encode the 
coin states. The transition from position \ak) to |afc±i) is achieved via two subsequent displacements, 
each followed by a pulse exchanging the the coin states (See figure |9]). (c) Probability distribution 
of the walker in position space after three steps, under the assumption that the position states are 
orthogonal. The black (white) filled boxes represent the contributions of the wave function related to 
|r) {\H)). The asymmetry between the position state probabilities P{ai, ip^) — [ij]^] (l® \ai) (ail) IV's) 
and P(a_i,'03) is due to the interferences indicated via the bold and dashed arrows in (a). 



according to 

T./. ,N ( cos {6/2) e^<^sm(^/2)\ ^ 

From the initial state l^/'o); the operator C with cj) being arbitrary, but equal for 
every application of C, leads to an asymmetric QW (Figured]). A symmetric QW can 
be realized with the coin operator for the first step being i? (71/2, 0) and for all following 
steps being R (7r/2, + 7r/2) (with arbitrary). In that case, the first coin toss can be 
interpreted as the initialization of the coin state such that all following coin tosses act 
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symmetrically on it. 

The shift operator 5* is defined as 

S = \T) (T| ® D{Aa) + \H) {H\ ® D(-A«), (7) 

with D{Aa) = exp(Aa-a^ — Aa*-a) being the displacement operator and a'^ ,a the 
corresponding raising and lowering operators. 

In contrast to a typical the position states \ak) are not orthogonal. The step 
size |Aq;| determines the overlap of the position states, {ak\ai) = exp(— (fc — /)^| Aap/2). 
If the state of the walker after steps is \iPn) = '^k=-N ('^fc 1^) I'^fc) + \H) l^fc)), the 
probability of finding the walker in position lai) is given by 



N 



Yl (^k-{oik\aL) 



k=-N 



N 



k=-N 



(8) 



Thus only if the step size |Aa| is large enough such that the overlap between different 
position states remains negligible, the above probability is given by the coefficients 
and c[ only. In figure [2] the probability distributions after 100 steps for QWs with 
different step sizes |Aa| are illustrated. We find that for |Aa| > 2, where the position 
states contain a negligible overlap of |(afc|Q;fc+i)p < e""^, the probability distribution 
shows the shape of an orthogonal QW [4j. For smaller values of |Aa| the probability 
distributions are smeared out due to the increased overlaps between the position states. 
As |Aq;| approaches zero, the probability distribution approaches a Gaussian shape. The 
mean distance of the walker from the origin, which is given by the standard deviation 



N 



(9) 



with 



Ef=-jv/(fe)--P(«fc,V'jv) 



ifik)) 



N 



(10) 



for any function f{k), grows slower with for smaller values of |Aq;|. This is because 
the lower the value of |Aa|, the less the shift operator S actually changes the state of 
the walker. However, for every realization of |Aa| the average distance of the walker 
from the origin scales linearly with the step number A^, i.e. ctat = f(AQ!)-A^, after a 
certain number of steps, as depicted in figure |2)3. The asymptotic scaling for the limit 
|Aa;| — )■ has not been investigated yet. 

In our experiment we set the step size to | Aa;| ^ 1. In this case the overlaps between 
the position states amount to |(afc|«fc+i)P = 1/^- The characteristic of the probability 



I In the sense that the position states are orthogonal, which has to our knowledge been assumed in 
the vast majority of publications concerning QWs, so far. 
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Figure 2. Characteristics of the probabihty distribution for asymmetric QWs of different step sizes 
|Aq;|. (a) Position probabilities -P(afc,V'ioo) after N = 100 steps. The probabihties are normalized to 
^f,P{ak,i^ioo)=^- For |Aa| > 2 the position states are approximately orthogonal and therefore the 
probability distribution shows the shape of an orthogonal QW 0]. In particular, the probabilities for 
positions \ak) with odd index k vanish. For |Aa| = 1 the probability distribution is smeared out and 
all position states are populated. However, the characteristic peaks around k = 70 are still prominent. 
For |Aq;| = 0.5 the probability distribution does not feature the main peak around k — 70 any more. 
For |Aq;| = 0.1 the probability distribution approaches a Gaussian shape, (b) Standard deviation ctn 
^ in dependence on N for QWs with different step sizes |Aa|. After a few initial steps the standard 
deviation scales linearly in N, i.e. crjv — w(|Aq!|)- N. For |Aa| = 0.1 already the initial state \ipo) of the 
walker is considerably spread out over the position states \ak) such that the standard deviation differs 
significantly from zero. The linear scaling becomes evident for N > 60. (c) Scaling factor ti(|AQ;|) of 
the standard deviation for different step sizes |Aq;|. For jAal > 2 the scaling factor is larger than 99% 
of the asymptotic value (w(|Aa| — > oo) = 0.457), which is the scaling of a QW with orthogonal position 
states. In our experiment we set the step size to |Aa| « 1, which results in a scaling factor of 89% of 
the asymptotic value. 



distribution after 100 steps (Figure [2^) is still close to that of |Aq;| > 2 and thus an 
orthogonal QW. 

Ideally, with the initial state \ipo) and = ([5]) three steps of the asymmetric QW 
lead to the state 

iV^s) = ^ |T) ® (-2 + |«3) - +^\H)^ (|a_3) + (H) 



The probabilities of finding the walker in the coin state \H) or |T) after three steps 
are given by Ph (V^s) = {H {\H) {H\ ® l) l^^) and Pt (^3) = {H {\T) (T| ® l) 1^3)- 
Their ratio 

Pnmi^l + e^^ .... 
amounts to ~ 1/3 for |Aa;| > 1. 
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3. Implementation of the QW 

3.1. System and definitions 

For the experimental implementation we confine a single ^^Mg"*" ion in a linear Paul trap 
|38] . The motional frequency related to the confinement in the axial direction of the 
trap is set to Uz = 2-7? ■ 2.13 MHz and in the radial directions to cj^; ~ cu^ ~ 27r • 5 MHz. 
We define two out of 12 electronic states of the hyperfine ground state manifolds [35] 
(Figure S]) 

\H)^\'S,/„F = 2,mF = 2), 
\T) = \'Sy2,F = 3,mF = 3) 

as the coin states. Further, we will use the state 

\A) = \^S^/2,F = 2,mF = -2) (14) 

in the detection procedure. To lift the degeneracy within each hyperfine manifold, we 
apply a magnetic field inducing a Zeeman shift with an energy separation related to 
Uzm ~ 27r- 3 MHz between neighbouring states. The energy difference (including the 
hyperfine splitting) between \H) and |T) amounts to a frequency of Ucoin = 27r- 1.77 GHz. 

Our realization of the QW consists of the application of a sequence of laser and 
radio-frequency (RF) pulses to (1) initialize the ion's electronic and motional state, (2) 
to implement the QW, and (3) to readout the final state via photon scattering. The 
experiments are repeated on the order of 1000 times for each set of parameters to obtain 
the required statistical relevance. A concise discussion of these tools in a generic context 
can be found in [39] and [iO] . 

3.2. State initialization 

At the beginning of each experiment the ion is prepared in the coin state |T) with a 
fidelity > 0.99 by optical pumping |4I], while the axial mode of motion is cooled close 
to the ground state by Doppler cooling (n ~ 10) [42] and subsequent sideband cooling 
(n < 0.03) [41j. The phase space of the axial mode of motion is used to encode the 
position of the walker. The radial modes are Doppler cooled {n ^ 4), which enables a 
sufficient decoupling from the axial degree of freedom. 

3.3. The coin operator 

We drive coherent transitions between the coin states by applying a radio-frequency 
field (RF) for a duration t with frequency Ucoin (Figures |3] and H]) [39j . This implements 
the operator R [6, 0) ([6]) with 6 = Q-t and Q = 27t- 100 kHz, the Rabi frequency of 
the transition. The phase (p for the first pulse of each experimental cycle is arbitrary. 
For every following pulse the phase <p is set identical with respect to the first pulse. 
The duration of one R (vr, 0)-pulse amounts to T^^ = 5 fis. The coherence time (drop of 
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Figure 3. Schematic of tlie operations required for the implementation of the QW. (a) Implementation 
of the coin operator R{d, cj)) ©• Left: The grey boxes represent a side view of two out of four electrodes 
of the Paul trap. The black dot between them depicts the trapped ion. We apply a radio frequency 
field (RF) via an antenna (below the grey boxes), driving coherent transitions between the coin states 
\T) and \H). The straight solid lines represent the phase fronts of the RF-field. Right: In our protocol 
of the QW we apply R{9 — ^, (/))-pulses {(f) = const.), rotating the state vector (Bloch representation) 
of \ipcoin) by 90° from |r) (|iJ)) to \H) + \T) {\H) - \T}). (b) State dependent optical dipole force for 
the implementation of the shift operator S ([7]). Left: We apply two laser beams, Rl and R2, with a 
frequency difference lui — uj2 = ujl — — perpendicular in polarization and beam direction. The 
effective wave vector is fc = fce^, pointing into the axial direction z. Right: This creates a walking 
standing wave and related state dependent ac-Stark shifts on the coin states \H) and \T), providing 
state dependent oscillating forces Ft (Fh) (solid sinusoidal lines) acting on the ion in the ^-direction 
with frequency lol- 



oscillation contrast below 50%) for the RF field exceeds one second [35j. The duration 
of a single experimental cycle (without initialization and detection) of the QW with 
three steps amounts to 150 /is (Figure [T0|) . Therefore the dephasing of the coin remains 
small. Spin-echo sequences |39], which are included in the QW pulse sequence (Figure 
[ro]) . further reduce the dephasing. 



3.4- The shift operator 

Ideally, we encode the positions of the QW into coherent motional states \ak) of the ion's 
axial harmonic motion. We manipulate the motion by implementing the shift operator 
S ([7]) via the application of a coin-state dependent optical dipole force (See figure |3b). 
In the following we describe this method and its limitations. 



3. 4-1- Experimental tools The initial motional state after sideband cooling is close 
to the ground state \n = 0) {n < 0.03). We apply a two-photon stimulated Raman 
transition between the coin states by applying two laser beams (i?l, R2) (Figures [3b 
andH]), at a detuning of A = 27r- 80 GHz from the P3/2 state manifold and a fixed phase 
relation [39]. The frequency difference between Rl and R2 amounts to 

= wi - W2 = ^2 - (5, (15) 

with 5 = 2TT- 100 kHz. The effective wave vector is k-^ — = ke^, pointing into the 
axial direction (Figure [3b). This allows for two-photon stimulated Raman transitions 
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Figure 4. Schematic of tlie relevant electronic and some of the lowest motional energy levels (not 
to scale) of ^^Mg'^ and the transitions used for the QW experiment. A magnetic field, pointing 
into the direction of the laser beam Rl, provides a Zeeman splitting of the hyperfine levels and the 
quantization axis of the system. The coin states are encoded in p5'i/2, ^"=3, mF=3) = |T) and 
\'^Si/2, F— "2,171^^2) ~ \H). Laser D, driving a closed cycling transition, in good approximation 
independent of the motional level \n), is used for optically pumping the electronic state into |T), 
for Doppler cooling of all motional modes 15] and readout of the internal state (See sect. 13. 5p |39) . 
The pair of laser beams BR and R2 is used to drive a two-photon stimulated Raman transition on 
the red sideband of the coin state transition for sideband cooling of the axiaconsequences of leaving 
the scope of the approximations originally used, such as the Lamb-Dicke regime. 1 motional direction 
[H] and on the blue sideband (BSB, shown here) for the readout of the motional state (See sect. 13.51) . 
The laser beams Rl and R2 drive a two-photon stimulated Raman transition providing the coin-state 
dependent optical dipole force (Figure [3)d). Additionally we apply a RF to drive coherent transitions 
between electronic states, independent of the motional state (Figure [3^). With the RF we implement 
the coin operation and the transition via several steps from \H) to \A) = \'^Si/2,F = — 2,tof = 2) 
required for the readout of the motional state (See sect. 13. Sp . 



|T) \n) ^ |T) \n + 1) and \H) \n) ^ \H) \n + 1) (Vn). In a simplified picture the two 
laser beams provide a walking standing wave causing a state dependent ac-Stark shift. 
This yields a coin-state dependent force {Ft, Fh), proportional to the spacial gradient 
of the walking wave and oscillating with frequency u^. The ratio of the forces acting 
on the coin states amounts to Fh/Ft ~ —2/3. The polarizations and intensities of the 
laser beams are adjusted such that the time-averaged ac-Stark shift for pulse durations 
T ^ 0.5 fis is negligible |13]. Thus the apphcation of the dipole force does not change 
the relation between the relative phase of the coin states and the phase of the RF, which 
implements the coin operator. The effective wavelength of the walking wave amounts 
to A ~ 200 nm. With the width of the axial ground-state wave function of zo ~ 10 nm 
this results in a Lamb-Dicke parameter of t] = zq- 27i / X = 0.31 [39] . 

3.4-2. Description of the dynamics We consider the following Hamiltonian describing 
a two-level system coupled to a harmonic oscillator and interacting with a classical light 
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field 



n ='H coin ~l~ T~^motion ~\~ 



motion ~r TLinteraction 



h 1 

^c^comO-z + f)uj^{a^a + -) + Mljj cos(fc- z - uit + 0o) (16) 



■ ~^^coin^z 



where fi^, = fi/j (|T) (T| — | with fi/j the coupling factor and the Pauli 

z-matrix. The dynamics have been investigated for many applications in QIP jHJ HS] 
(in the LDA, see below), for the simulation of nonlinear optics [16] and in the context of 
mesoscopic entanglement [331 |3l]. In the interaction picture, with the free Hamiltonain 
being 1-icoin + 'Hmotion, the interaction Hamiltonian can be written as 

. oo oo 

ni{t)= -Qd ®Y.Y1 \^)(^\ e''^^"^"^ ^\n){n\ 

m=On=0 (17) 
^ f ^i{{m-n)ui2-uiL)t+i<j)0 _j_ ^_^ym-n\ ^i{{m-n)LUz+LUL)t~i<f>o 



We apply the dipole force with a small detuning of 6 = — ujl = 2n- 100 kHz, such 
that the terms corresponding to first-sideband transitions, \n) ^ |n + 1), rotate slowest 
and thus dominate. However, we also consider contributions up to the third sideband, 
\m — n\ = 3, an approximation we refer to as 3SB (Figures El IHl and fT2l) . 

When considering only the slowest rotating terms, i.e. applying the usual rotating- 
wave approximation (RWA), the interaction Hamiltonian is reduced to: 

- oo 

nf^^it) =-Qo ® + 1| e'''("+"') \n) e'('^*+'^») \n + 1) {n\ 

2 „=o^ (18) 

- {n\ e*^''(»+'^^) \n + 1) e-^('^*+'^«) \n) {n + l\). 



For states with k\f(z^ = rj ^^/{{a+~aJy) ^ 1, we can approximate 
(n + l|e*''("+"^) \n) ^ iri\^Ti^\^ [ID]. That is, the potential providing the dipole force 
changes linearly over the extension of the wave function. We can then simplify the 
interaction Hamiltonian to 

nf'it) = ® zr/(ate*(^*+^") - ae~^('^*+<^«)). (19) 

This is the Lamb-Dicke approximation (LDA) [39] or the linear approximation, respec- 
tively. In the following, we set (po = and h = 1. 

Within the LDA, the time evolution operator is given by [17] 

U{t) = \T) (T| ® D (a(t)) ■ e'*("(*)'*) + \H) {H\ ® D (-^« " e'*(^i"(*)'*) , (20) 
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Figure 5. Normalized Fock-state transition rates |fin+i.n|/|i^i.o| — K"- + 1| exp(i77(a + a''')) |ri)| p8|) 
for a Lamb-Dicke parameter rj = 0.31 (experimental value) and rj = 0.001 (deeply within the LDR). 
As long as |f2„+i^„|/|rii^o| ~ '>lV^^ the force is constant over the ion's motional extension. The time 
evolution is then described by a displacement operator ((20)) . preserving the shape of coherent states. As 
the motional amplitude |a| = \/Jnj increases, the force starts to remarkably change over the oscillating 
range of the wave function leading to motional squeezing. The Fock state at which the matrix element is 
maximal, i.e. n = 8 = gi for rj = 0.31, can be considered as the threshold above which severe motional 
squeezing of coherent states starts. The motional amplitude does not increase by the application of a 
continuous detuned force (See figure IH). Motional states up to 52 for 77 = 0.31 can be created - with 
an increasing amount of squeezing - by applying the dipole force resonantly (Figure [7|) or by step- wise 
excitation (Figure [8|). 




which is a displacement operator D {a{t)) with a phase factor, where the phase amounts 
to 

da{T) 



$(a(t),t) = Im (^J^ dTa*{T) 



dr 



(21) 



The factor 2/3 in the displacement operator for the Heads-part results from the 
difference of the state-dependent dipole force, Fh/Ft = —2/3 (See figure [3]). The 
complex parameter appearing in the displacement operator and in the phase amounts 
to 



aW = ^-(e«* = -.!§^.(e'«-l) (22) 



and corresponds to a circular trajectory in a co-rotating phase space, given by the 
interaction picture as 



Re(a(t))\ ^ / cos(w,t) sin(a;,t)\ ( (z) {t) 
Im(a(t)) I I - sm{u,t) cos{u,t) \ Zq {p) (t) 



(23) 



with (z) (t) and (p) (t) the expectation values of position and momentum. 

For each coin state, the motional wave function is coherently displaced along a 
circular trajectory in the co-rotating phase space (Figure E]). The circular shape of the 
trajectory is caused by the dipole force being applied with a detuning S relative to the 
oscillator frequency. Thus the relative phase between dipole force and oscillation of the 
ion evolves in time as (pnit) = ^- 1. After a duration of = 7r/5 of driving the motional 
state and increasing its amplitude, the relative phase amounts to 0z)(T^) = vr and 
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Figure 6. Numerical simulation of the ion's trajectory in co-rotating phase space ((23)) driven by the 
dipole force for the three relevant approximations 3SB P7|) , RWA (IT51) and LDA . The initial state 
is at the origin (|ao — 0)). The thin concentric lines represent contours of its Wigner function W (at 
VF> = 0.6 and VF< = 0.3). The bold dotted line represents gi (Figure [5]). The thin circular trajectory 
with dots represents the result of the simulation within the LDA. The dots on the trajectory depict 
the positions after t = 0,0.5, ...10 jis. The final state, reached after = 2tt/6 — 10 jis, equals the 
initial one, up to a phase factor. The bold trajectory represents the result within the RWA, taking 
nonlinearities of the dipole force into account (Figured]). The dots on the trajectory again depict the 
position at the times t — 0,0.5, ...10 ps. Starting from the origin, the trajectory is identical to the 
one within the LDA. In the vicinity of gi the trajectories start to deviate. The acceleration of the ion 
ceases at a certain amplitude, the state gets squeezed and then returns to the origin after a duration 
shorter than T2Tr- The spiraling trajectory, which follows the one within the RWA, represents the 
results within 3SB. Here terms of higher frequencies in the Hamiltonian are taken into account. The 
final Wigner function is almost identical to the one in the RWA and therefore not shown. Parameters: 
fir, = 27r • 1.2 MHz, wl = 27r • 2.03 MHz,w^ = 27r • 2.13 MHz, r] = 0.31, t G [0,10] jis. 

therefore the dipole force starts to decelerate the oscillation of the ion. After a duration 
of T27r = 271/6 the coherent state returns to its initial location in the co-rotating phase 
space. The total acquired phase of the motional state, which equals the enclosed area 
of the trajectory [H], amounts to 

for |T) and $h = (4/9)- $t for \H). 

The nonlinearity of the potential causing spatial variations of the dipole force can 
be described by the absolute values of the transition matrix elements, i.e. fln+i,n = 
(n + 1| exp(^ir]{a + a^)^ \n) (Figure [5]). For small Fock state numbers n they remain close 
to the approximative results within the LDA (where the potential is linear in z). Close 
to n = 8 = Qi (for rj = 0.31) they start to significantly deviate from that approximation. 




(24) 
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In particular, the transition matrix elements feature a maximum value at gi. The 
excitation of motion via the dipole force ceases at n = 37 = g2, due to 1^38 37! ^ 0. 
It therefore represents an upper bound for the motional excitation with a dipole force 
applied close to resonance. (Applying the dipole force with a frequency ul = '2-oJz 
would allow populating higher Fock states, but since the overall time evolution is then 
described by a squeezing operator, it cannot be used for the implementation of a QW 
based on coherent displacements, following reference [28].) 

Figure [6] presents the results of our numerical simulations of the time evolution, 
comparing the three different approximations 3SB, RWA and LDA. Starting in the 
motional ground state, the trajectory in the co-rotating phase space first follows the 
circular evolution, as long as the amplitude remains small, i.e. (n) < gi. This regime 
can be well described by the LDA (with the driving potential being linear in z). As the 
amplitude of the motional state approaches gi, the driving potential becomes sufficiently 
nonlinear to severely affect the subsequent evolution. The amount of displacement per 
time interval is substantially reduced, as the transition rates |n„+i,„| decrease for n > gi, 
such that the trajectory remains in the vincinity of gi. At this point motional squeezing 
occurs. The probability distribution of the wave function in the Fock state basis becomes 
narrower than Poissonian, which results in a squeezed shape of the corresponding Wigner 
function. The relative phase between the dipole force and the oscillation of the ion 
changes faster than in the linear case. Thus the squeezed wavefunction reaches the 
origin of the phase space after a time significantly shorter than 2tt/6. The dependence 
of the return time and the amount of squeezing on the maximal motional amplitude 
severely affect an implementation of a QW with position states outside the LDR, 
following the scheme described in reference p8j. However, in section [6] we propose 
an alternative protocol for the implementation of the shift operator that circumvents 
this restriction. The faster rotating terms, which are taken into account in the 3SB 
approximation, cause additional modulations of the trajectory with low amplitudes and 
high frequencies {2uz + S) and {Su^ + S), respectively. 

Further motional excitation (up to (72) can be realized either by applying the dipole 
force resonantly (Figure [7]) ^46], or by the repeated off-resonant application of a weak 
dipole force with duration n/S and constant time delays between the pulses (Figure [8]). 
However, in either way severe motional squeezing occurs as the amplitude approaches 
and becomes larger than gi. 

3.4-3. Implementation of the shift operator We implement several shifts into a certain 
direction in the co-rotating phase space by a synchronized application of dipole force 
pulses. Switching the dipole force is realized by acousto-optical modulators, refracting 
the laser beams Rl and R2 into the Paul trap [39] . As the laser is continuously on during 
the whole experiment, the phase relation between the dipole force and the motion of the 
ion continuously evolves in time as (J) nit) = S-t, even when the optical dipole force is 
switched off. Therefore, after the first pulse of the dipole force, the relative phase (f>D{t) 
and thus the direction of the displacement caused by the following dipole force pulse 
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Figure 7. Numerical simulation (3SB) of a resonant (5 = 0) excitation from the motional ground 
state, trajectory (grey horizontal line) and contours of the Wigner function W (black: W = 0.1, grey: 
W — —0.1) for increasing durations t = {0, 1.5, 2, 4, 6, 8} jis. Starting from the ground state, the dipole 
force continuously displaces the ion in good approximation to coherent states along the real axis up to 
gi (dotted circle). Further excitation comes along with severe motional squeezing. Amplitudes higher 
than g2 are not considerably populated, due to the almost vanishing matrix element 1^37. as at 32 (dash- 
dotted circle) (Figure [S]). As a consequence, at 52 the wave function gets reflected in the sense that 
the amplitude gets decreased again. The interference of the accelerated and the decelerated (reflected) 
part of the wave function is resembled in the Wigner function, which shows concentric lines of positive 
and negative value in the populated area of the rotating phase space. Parameters: = 27r • 2.0 MHz, 
= a;^ = 27r • 2.0 MHz, 7/ = 0.3. 



depends on the intermitted delay. We apply pulses of the duration T5 ^ 7r/(5 = 5 ]js 
and mutual delays that concatenate the displacements along a line in the co-rotating 
phase space (Figure [8]). 

However, the motional states corresponding to \H) and \T) acquire different 
phase factors $r (12^ during each displacement. To compensate these coin-state 
dependent phases, we implement the shift operation of the QW as a combined pulse, 
which consists of two dipole force pulses, each followed by an i?(7r, 0) -pulse (Figure [9]). 
In this scheme each coin state acquires the sum of the phase factors, $7 + ^h-, which 
therefore turns into a global phase factor not affecting the QW. A schematic of the 
overall pulse sequence for the QW is depicted in figure (ITO!) . 

For the implementation of the QW the following is crucial. After a step fulfilling 
the operation |T) |5fc) — )■ |T) \oik+i), the subsequent coin toss and shift operation have 
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Figure 8. Numerical simulation (3SB) of the ion's trajectory in the co-rotating phase space during a 
stepwise excitation, (a) The first dipole force pulse (inlay) displaces the ion from the ground state to 
|Si.) (which is not the position state |ai), compare figure[9]) along trajectory ti. During the pulse the 
relative phase between dipole force and oscillator motion increases to (j)D = tt. The phase difference 
further increases after the dipole force is switched off (Sect. 13.4.31) . After the duration t = tt/^, it 
amounts to 0£) = 27r. Applying now a second dipole force pulse (inlay) displaces the ion to |52*) 
along trajectory T2. (b) Stepwise excitation as described above, with eight steps, passing the threshold 
g\. As the state passes gi, the trajectory during each dipole force pulse follows the opposite direction 
of rotation (counter-clockwise). In this regime (in particular close to g-z) the time evolution can be 
approximated by a displacement operator as in (|20|) parametrized with a(--t) 46 . The contours of 
the final Wigner function show the strong motional squeezing. In addition, due to the dependency of 
the return time on the amplitude (Cf. Figure [6]) the trajectory is on average not oriented along the 
horizontal axis. This is also the case in the protocol of the three-step QW f Figure [T0| . but does not 
affect its performance. Parameters: fi/, = 27r • 0.4 MHz, = 27r • 2.0 MHz, 8 = 2-k- lOOkHz, -q = 0.3. 



to ensure the operation \H) |5fc_|_i) — )■ \H) \ak), i-e. the motional state \ak+i) must be 
transferee! back to the previous state \ak)- This must be fulfilled for all k simultaneously 
(See figured]). In order to reach the state |5fc), the duration To and detuning S of the 
dipole force have to be adjusted properly (In the LDR to To = tt/S), implementing 
semi-circular trajectories in the co-rotating phase space. But since the return time is 
reduced outside the LDR (Figure [6]), the shift operation implements the transition from 
|5fc+i) to some other state IS^j) 7^ |5fc). The reduced overlap (SfcjISfc) < 1 leads to 
reduced interference during the succeeding coin toss (Figured]). 

In principle the shift operator can alternatively be implemented by the dipole force 
on resonance {6 = 0). However, in that case small variations of and Uz have a 
much stronger infiuence than in the detuned case. This can be seen by comparing the 
difference in displacement, i.e. \a{6,t) — a{6 + e,t)\, for a fixed duration t and a fixed 
difference in the detuning, e, using ( 122]) for different values of S, where a{6,t) fulfills 
equation ( 122]) for the detuning 6. 

3.5. State readout 

The state after three steps of the QW is 

= (|T) \Mt) + \H) \Mh)) , (25) 




Figure 9. Numerical simulation (3SB) of the trajectories related to the ion's coin states |T) (bold 
trajectory) and \H) (thin trajectory) during step 1 of the QW fFigure [T0| . After the coin operation 
the ion's state is lip) = (|r) + \H)) \ao) (where \ak) denotes the possibly slightly squeezed version of 
|q;a;)). The shift operation S is implemented by a combined pulse. The first dipole force pulse displaces 
the motional state related to |T) {\H)) to \ai*) (|S_i»)). As the forces are of different amplitudes, 
the two different trajectories lead to different phase factors, related to $t, ([Ml)- Subsequently the 
coin states are exchanged via an R{tt, 0)-pulse without affecting the motional states and after a specific 
waiting duration (Figure I10[) a second dipole force pulse is applied, displacing the motional states to 
|5_i) (|5i))- A second i?(7r, 0)-pulse exchanges the coin states again, such that the resulting state of 
the ion is \^pi) = exp(i($T + ^h)) - {\T) |ai) + \H) |a-i)), accumulating only a global phase during the 
shift operation. The combined-pulse scheme further provides equal step distances in both directions. 
Parameters: no = 2tt ■ 0.24 MHz, W;, = 27r • 2.13 MHz, 5 = 2tt ■ lOOkHz, r/ = 0.31. 



with \Mt) = Ylk=-s^k I'^k) and \Mh) = Ylt=-3'^k I'^fc) (^f- dH]))- The basics for the 
readout are state-of-the-art techniques in quantum information processing with trapped 

ions llHlEg]. 

We readout the coin state by driving the cychng transition |T) |?7,) — )■ 
I2P3/2, -F = 4, mi? = 4) |?T,) (Figure Hj) for a duration of 20 ]is and detect scattered 
photons with a photomultipher. This transition is in good approximation independent 
of the motional state. The average number of detected photons is proportional to 
the probability Pxiips), related to the coin state |T) [39]. The probability Pni'ips) is 
accessible via the application of a i? (tt, 0)-pulse before the coin-state detection. 

To characterize the motional states of {ip^) we determine the position-state 
probabilities jcl'p and |c^p. To analyze |T) \Mt), we isolate the other part, \H) \Mh), 
by transfering \H) \Mh) \A) \Mh) using appropriate RF pulses (Figure H]). The part 
1^4) \Mh) of the ion's state is not affected by the subsequent operations. We then apply a 
two-photon stimulated Raman transition |T) \n) -H- \H) \n + 1) (Vra) using the lasers BR 
and R2 (Figure H]), with a frequency difference of ojl = ojcoin + (BSB) for a variable 
duration Ibsb |39]. The corresponding Rabi frequency for each n is proportional to 
VLn+i,n (Figure E]). Finally, we apply the coin-state detection, as described above. The 
average number of detected photons is proportional to the probability [39] 

PxitBSB) = i ( 1 + X]|a^|2. cos(r]„+i,„, ■ W)e-^ *^^^ ) , (26) 

V n=0 / 

{n\MT) the coefficients of \Mj) in the Fock state basis. The damping 
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Figure 10. Schematic of the pulse sequence for the implementation of three steps of an asymmetric 
QW. The first line (top) illustrates the overall pulse sequence for one experimental cycle, in chronological 
order from left to right. We initialize the ion in the state \tpQ) — \T) \ao). Then we apply the pulse 
sequence for the QW, creating the state IV'3) (|25p . and finally a pulse sequence to readout the coin or 
motional state (Sect. 13. 5|) . The pulse sequence for the QW consists of three subsequent applications 
of coin C (Figure [3]) and shift S (Figure E]) operations, implementing the three steps of the QW as 
illustrated in the second line. The third line depicts the pulses implementing one coin and one shift 
operation. The phase 4> (not depicted) is equal for every RF pulse (Sect. 13. 3|) . The sequence of all RF 
pulses incorporates a spin-echo scheme, reducing dephasing of the coin states. A symmetric QW can 
be implemented with the same scheme, but with the phase of the initial RF pulse differing by 7r/2 from 
the (equal) phases of all other RF pulses. The timing of the pulses, in particular for the dipole force, 
is given by the parameter Tjj , as illustrated in the bottom. It must be chosen such that the overlap of 
the interfering motional states (Figure [1]) is maximal. The theoretical optimal value within the LDR 
would be Td — tt/6. We determine the optimal value of To, denoted as sa 7r/<5, experimentally to 

account for the nonlinearity of the dipole force as well as for experimental imperfections (Figure [T2| . 



factor 7 accounts for decohering effects, mainly of the BSB-operation |39]. A discrete 
Fourier transform of the function PritBSB) allows to access the Fock state probabilities 

We numerically simulate the QW and generate a corresponding function Pj'itBSB) 
(where the index ^ denotes the result of the simulation), and optimize the parameters of 
the simulation in order to fit PxitBSB) onto the experimental data, P^itBSB) (indicated 
by the index ^). The state j-^f), generated by the simulation, resp. the Fock-state 
probabilities of its motional parts, are fed into the algorithm for the discrete Fourier 
transform of P^itBSB) to improve the convergence of the results. The Fock state 
probabilities Pn (5^) = |(n|5fc)p of the position states |5fc) are determined separately 
(Figure [TT]) . using the same method. 

To identify the position state probabilities |c^p from the Fock state probabilities of 
I Mr), in particular to distinguish between the coefficients related top„ (5^) andp„ 
we additionally apply the motional-state readout to the shifted states l?/'^) = S \ipz) and 
I'tp^) = lips) (See ([7]) and figurelH]), where is (up to a global phase) implemented 
by a proper timing of the corresponding dipole force pulses. 
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Figure 11. Fock state probabilities p„ (af) = of the simulated (3SB) position states 

resembling the probabilitites p„ (5f^ ) of the experimental position states | af ) . Experimentally, we apply 
the shift operation S (Figure IH]) k times (without coin operations) to the initialized ion and readout 
the motional state (Sect. 13. 5p . The inlay shows the coin state probability P^{t) (black squares) from 
the motional state readout of jT) jaf), as an example. Each experimental data point represents the 
average of 3000 realizations. A simulation (3SB) of this procedure is optimized such that the values 
P^{t) (dotted line, cf. (|26)) ) correspond to the experimental result for every position state |af), as 
shown for |af ) in the inlay. The related expectation values of the position states in the Fock basis 
(main figure) are (n )o = 0, (n)l = 1.33, {n)l = 4.71, {n)l = 9.08, (n)^ = 13.50. Their mutual overlaps 
amount to |(af |af_^^)|2 « 0.24. 



4. Experimental procedure 

J^.l. Calibration of the step size |Aa| 

Starting from the initial state we apply to 4 shift operations S with a dipole force 
duration of ~ |, exciting the motion of the ion to one of the position states |5fc). 
Then we readout the motional state to determine its Fock state probabilities pn iptk) 
and the expectation of the number operator (n) = |5fcp. Small deviations of the dipole 
force duration (See next subsection) do not influence the probabilities significantly. We 
adjust the amplitude of the dipole force by adjusting the corresponding laser beam 
intensities to approximately meet the conditions |Aa| > 1 and ISsp = ifit^\n\a^) < 9 
(three well distinguishable position states within or close to the LDR). The Fock state 
expectation values of the position states amount to (n)^ = (5o| n |5o) = 0, {n)^ = 1.33, 
(72)2 = 4.71, (72)3 = 9.08, {n)^ = 13.50. The outer position states |5±3) and |5±4) are 
not within the LDR. The step sizes therefore differ, reaching from |5i| — |5o| = 1.15 
to I54I — ISsI = 0.66. However, due to the motional squeezing the overlaps of all 
neighbouring states amount to |(5fc|5fc+i)p ~ 0.24 < 1/e, which corresponds to the 
overlap of coherent states with a step size of |Aa| > 1 as assumed in the theory (See 
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Figure 12. Coin state probabilities in dependence of the relative dipole force pulse duration {T]:)/{j)) 
after the application of the QW pulse sequence (Figure lTOl) to the ion in state | V'o) : and the corresponding 
numerical simulations (solid line: 3SB, dashed line: RWA). Each data point represents the average of 
1500 realizations. The maximum ratio Ph / Pt = 1/3 indicates the asymmetry due to interference (See 
figure [T]). The corresponding dipole force duration, denoted as T^^, is the optimal value to perform 
the QW. The precise value of depends on the detuning 5 which is prone to slow drifts of the 

conditions of the experimental setup (on a time scale of a few hours, much longer than an experiment) 
and can be estimated to the required precision (for each experiment) by this method. For other values 
of To (i.e. 2% longer or shorter) the effect of the interference vanishes, since the shift operations of 
the pulse sequence do not lead to mutual overlaps of the related parts of the wave function. In this 
experiment the QW pulse sequence (Figure [T0| contains waiting durations of ATe, instead of 2Td, which 
increases the sensitivity of the interference to Tu and therefore allows for a more precise estimation 
of the optimal dipole force duration T^^. The 3SB simulation (solid line), in contrast to the RWA 
simulation (dashed line), contains a high-frequency modulation of the coin state probabilities due to 
additionally modulated overlap of the interfering states caused by the spiralling trajectories during the 
displacements (Figure [6]). This is however not yet resolved in the experimental data. 

section [2]). With 5 = 2ti- 100 kHz and rj = 0.31 the corresponding couphng strength of 
the dipole force amounts to Qd = 27r-0.24 MHz. The amphtude of the corresponding 
dipole force amounts to Ft = [hrjVtjj) / {2zo) = 2.54- 10~^^ N inside the LDR. 

To estimate the amount of motional squeezing (not affecting the fidelity of our 
results on the QW), we compute coherent states \ak) according to equation ([3]) with 
ak = 5fc. We then compute their overlaps = |(5fc|afc)p, which amount to Fq = 1.00, 
Fi = 1.00, Fa = 0.97, F3 = 0.90, F4 = 0.78. 

4-2. Calibration of the dipole force duration 

The duration To of the dipole force pulses (and the related mutual delays, see figure [TOll 
is a very sensitive parameter for the implementation of the QW. For a given detuning S, 
Td determines the relative direction of subsequent shifts in the co-rotating phase space. 
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Figure 13. Numerical simulation (3SB) as depicted in figure [T^ for extended values of Tq. Next to 
the largest difference between Pt (^ipi^ and Ph (^'Af ^ at To « 5.0 jis ~ tt/S, similar splittings occur at 
Tjj « 4.6 jis and « 5.4 jis. These durations are such that subsequent shift operations displace into 
opposite directions in the co-rotating phase space. That is, at odd step numbers of the QW the shift 
operation acts as \T) \ak) — > |r) and at even steps as |T) \ak) — > |r) \ak-i), and analogously for 

\H) |afe). This creates again a QW in which the shift directions of |T) and \H) are exchanged at each 
step. 

By altering we can control and maximize the overlap of the interfering parts of the 
wave function in the QW (Figure [1]). We repeat the QW pulse sequence with increasing 
values of To (Figure [TO]) and acquire the coin-state probabilities Pt and Ph via the 
coin-state detection (Figure [T2|) . The ratio Pt/Pr is an indicator for the amount of 
interference. If no interference occurs, it amounts to Pt/Ph = 1- We maximize the 
ratio by iterating to the optimal dipole force duration, denoted as ^ ti/5. The 

nonlinearity of the dipole force, in particular the reduced return time discussed in sect. 
I3.4.2[ leads to a deviation of from the duration 7r/5, optimal within the LDA 

only. Additionally, this method is suitable to implicitly determine the detuning 5 to the 
required precision. 

As illustrated in figure [121 the maximum ratio of the measured probabilities 
amounts to Pt/Pr ~ 3, with the related dipole force duration being defined as . 
Averaged over 60000 measurements the coin state probabilities at this point amount 
to Pt = 0.741 ± 0.002 and Pr = 0.259 ± 0.001, which is close to the theoretical 
predictions 0.75 and 0.25 (Figure H]). At slightly different values of the dipole force 
duration, Td = T§^- (1 ± 0.02) the coin state probabilities are approximately equal 
{Pt/Pr ~ 1), indicating that the overlaps and hence the interference of different parts 
of the wave function vanish. 

The results of a numerical simulation of this procedure within 3SB are in good 
agreement with the experimental data (Figure [T2|) . Additionally, the simulation shows 
similar splittings of the coin state probabilities at the dipole force durations Td = 4.6 fis 
and Td = 5Afis, for 6 = 2tt- 100 kHz (Figure [T3l) . These are QWs in which the step 
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Figure 14. Fock state probabilities Pn{ip) = \{T\ {n\ip)\'^ and = \{H\ {n\ip)\'^ of the simulated 
QW state I'^f ) (top), the shifted states ) , {-tp^" ) (center, bottom) and of the position states jaf ) 
(dotted lines), as in figure [11] All simulations are computed with equal parameters (T^^ = 4.990 jis 
and the 3SB approximation). Each simulated state is in agreement with the corresponding experimental 
data, as exemplarily depicted in figure [TT] (inlay). Top: lip^) contains a superposition of \ai) and |a_i). 
This causes an interference (not related to the interference of the QW) in Fock space, resulting in high 
probabilities for even {n = 0,2,4) and low probabilities for odd Fock states. With the states li^^'^) 
(center) and \ip^' ) (bottom), where the position occupations are shifted by one position, it is possible 
to distinguish the probabilities corresponding to \ak) and |a_fe). 



directions of \T) and \H) are exchanged at each step. 

The simulation also shows a high frequent modulation of the coin state probabilities 
in dependence of To. This is due to the fast modulations of the trajectories in the 
rotating phase space, mainly caused by |n + 2) -H- |n)-contributions of the dipole force 
(See figure [6]). 



Experimental Simulation and Limitations of QWs with Trapped Ions 



23 




position state 

Figure 15. Position probabilities P{a^,ijjf) after three steps of the asymmetric QW. These are 
acquired by a discrete Fourier analysis of the data from the motional state detection (|26p of the states 
Tpf, V'a"' J '03~' (Cf. figure [13]). The error bars represent the errors of the discrete Fourier analysis. 
The probabilities are in agreement with the theoretical values of an asymmetric QW (See figure [T]). In 
particular the difference of the probabilities P{af,^f) and P{a^i,tp^) indicates the high fidelity of 
the implementation of the QW. The probabilities at the positions \a^2)' I'^o) ''^^'^ I'^f) rema-in nonzero 
due to the overlaps with the neighbouring position states. The probabilities due to the overlaps have 
been subtracted, using the probabilities at \a^i) and \a^^) as a reference. Therefore the remaining 
nonzero values are an indicator for the error of the implementation and readout of the QW. 

5. Experimental results and conclusion 

We implement the QW pulse sequence (Figure [10]) with the optimized dipole force 
duration (Figure [I2]) and apply the motional state readout (Sect. 13. 5p . The 

resulting Fock state probabilities from the corresponding simulation are illustrated in 
figure [HI To distinguish the position states |5fc) and |5_fc), which have the same Fock 
state probabilities, we also apply the motional state readout to the states {ipt) and l^p^), 
where after the QW pulse sequence an additional shift operation towards higher (lower) 
position states has been applied. The position state probabilities corresponding to the 
experimental data are illustrated in figure [151 

The experiment demonstrates the feasibility of implementing a QW with a trapped 
ion. Although the number of steps is small in our experiment, the trapped- ion system 
clearly reveals its strengths in the high fidelity of the results. The severe limitation on 
the number of steps for the implementation of the shift operator via the optical dipole 
force is due to the Lamb-Dicke parameter rj, since shifts in terms of the displacement 
operator are only possible within the LDR (Figures |5] and |6]). Our comparatively large 
Lamb-Dicke parameter rj = 0.31 allows for three well distinguishable (|Aa| = 1) steps. 
For a Lamb-Dicke parameter of rj = 0.1 the limit of the LDR is gi = 85 with a 
corresponding maximal position state \amax ~ 9) within the LDR. This would allow for 
a QW with 9 steps of similar fidelity, using our scheme. Equivalently, a Lamb-Dicke 
parameter of rj = 0.06 would allow for 15 steps. Additionally, since in such a setup 
the step size is very small compared to the size of the LDR, the threshold gi may be 
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less of a limitation, as it can be overcome via small steps (Cf. figure [H])- A QW with 
up to 23 steps, implemented with a dipole force on resonance {6 = 0), has recently 
been demonstrated [32]. As described in section Ej a QW with effectively orthogonal 
position states requires a step size of |Aa| > 2. This reduces the number of steps within 
the LDR for any given setup. Extending the number of possible steps substantially by 
further reducing the Lamb-Dicke factor is a difficult task. The trap frequency has 
to be increased and the mutual angles of the laser beams providing the dipole force 
must be reduced (Cf. figure [3b). However, a small Lamb-Dicke factor rj yields a weak 
coupling of the light field to the motional degree of freedom. This weak coupling must 
be compensated by an increased intensity of the laser beams. This in turn results in an 
increased rate of spontaneous emission [13] from the off-resonantly excited P-P3/2, F = 4) 
states (Figure Hj) and thus in a reduced coherence time for the QW. 

6. Outlook 

6.1. Implementation on the shift operator using photon kicks 

In the following we propose the implementation of the shift operator with photon kicks 
[361 [37] , which is substantially less dependent on the motional state and allows for the 
implementation of QWs with many steps. The principle of a photon kick is to apply a 
vr-pulse on the coin states which is sufficiently short such that the free harmonic motion 
of the ion during the pulse itself is negligible. It has been shown that the change of the 
momentum of the ion during such a pulse can be described by a displacement operator, 
allowing us to propose its application as a building block for the shift operator of a 
QW. In the original protocol |36j, which has been realized recently [19], however, the 
infiuence of the motional state on the performance of the photon kicks has not been 
considered, since the amplitudes of the motional states were assumed to remain small. 
For the implementation of a QW with many steps, we have to consider (coherent) 
motional states with very large amplitude and, thus, have to re-assess the validity of 
the above-mentioned approximation. We find that, for a given fidelity, the upper bound 
for the pulse duration scales inversely with the motional amplitude, and additionally, 
for coherent motional states, depends on the phase of their harmonic oscillation at the 
moment when the pulse is applied. In the following we derive an analytic bound for 
general states and present the results of a numerical study for coherent motional states. 
With the latter we show that QWs with up to 100 steps for a step size of |Aa;| = 2 
should be possible with state-of-the-art technology. 

Referring to [36], we start our analysis with the Hamiltonian 




This Hamiltonian can be implemented in various ways, e.g. via direct dipole 
coupling, two-photon stimulated Raman transitions or stimulated Raman adiabatic 





(27) 
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passage Each implementation imposes different constraints on pulse duration, laser 
intensities, etc. In the following we will focus on the implementation with a two-photon 
stimulated Raman transition and consider the energy levels of ^^Mg^. 

In this configuration, two laser beams {Ra,Rb) resonantly drive two-photon 
transitions between the coin states via a virtual state detuned from the P3/2 manifold by 
A/j. Each laser beam drives only one of the two Raman branches, due to their different 
polarizations. In a RWA, terms varying at optical frequencies are neglected. This is 
valid in our case for pulse durations well above 1/10^^^ Hz = 1 fs. Finally, an adiabatic 
elimination of the P3/2 states requires |f2/Aj:j| ^ 1. The pulse duration Tp in our case 
must therefore be sufficiently longer than 5 ps for ^ 2tt- 10^^ Hz and TpQ = tt (see 
below). The effective wave vector of the two-photon transition is k = ka — kf,. 

Hamiltonian ( l27l) implements the desired displacement operator for a pulse duration 
of Tp = tt/Q, if we neglect the perturbation "Hi. The time evolution operator then reads 



which is obtained by expanding the exponential function, splitting the series into odd and 
even parts and using the properties of the Pauli matrices and displacement operators. 

The shift operator itself, implementing the desired step size (i.e. |Aa;| = 2, see 
figure |2]), can be realized by the subsequent application of 2/?7 kicks in such a way that 
the displacements Diirf) of several 7r-pulses add up to ^(Aa = iri-2/ri). This can be 
achieved by changing the direction of the effective wave vector by 180° for each photon 
kick. In practice one can either switch between two Raman beam configurations with 
opposite effective wave vectors, or implement every second vr-pulse by a RF transition 
for which the momentum transfer is negligible. Notably, with this protocol the step 
sizes for both directions of the QW are equal, in contrast to the method of optical 
dipole forces used in our current experiment. 

In the following we derive a conservative estimate for the deviation from a coherent- 
state displacement induced by T-ii. The total time evolution is 



UoiTp) =e 



-moTp 




- i sin (^-^J ■ (0-+ (g) D {ii]) + a_®D {-irj)) 
-i{a+®D {ir]) + a^®D {-ir])) , 



(28) 



U{t) = e 



-mt 



Uo{tyv{t) 



(29) 



with V{t) = e'^ot^-iHt_ y(^^^ ^g^^ differentiated 

V{t) = -ie'^'''nie-''^''-V{t), 



(30) 



leading to an equation wich is formally solved by the integral equation 




(31) 
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using V{0) = I. Now we define e as the distance between the evolved state according 
to the full Hamiltonian and the desired evolved state according to T-Lq: 

e^ll {U{t)-Uo{t))\^)\\ 
= II {V{t) - I) \^) II ^32^ 

= II / ds e'^^'Hie-'^'^'Vis) |^) ||. 
Jo 

Approximating the last expression by the largest term of the first order Dyson series, 
and considering a motional state = \H) \a), gives the following error estimate for a 
pulse with duration Tp [50] 

ds cOzl ® a^a \\ = TpUJz\a\'^- (33) 







Thus, for an initial state \H) \a) (the coin state can be chosen arbitrarily), the pulse 
duration Tp necessary to implement the displacement operator with an error smaller 
than e must fulfill 

Tp < (34) 
(^z- |ar 

The scaling with |a|~^ is, however, a rather rough estimate. This is shown by 
a numerical simulation of this process, in particular considering the application of 
photon kicks to (superpositions of) coherent motional states. We compute the fidelity 
/ = \{H\ {a\UQ{Tp)U{Tp) \H) |q;)P with the initial state \H) \a), a pulse duration Tp, 
and Q = ir/Tp, where the time evolution is implemented using a Runge-Kutta method. 
The results show that the fidelity strongly depends on the phase of ion oscillation at the 
moment of the photon kick. 

Demanding a fidelity of / > 0.99 and for imaginary a, i.e. at the moment of the 
photon kick the ion is in the center of the harmonic potential and thus fastest, for our 
experimental parameters we find [§| 

^/Lo.99(l«l) = exp (-17.55 - 0.631n(|a|) - 0.05 (In (|a|))^) . (35) 

For |q;| = 200, an amplitude reached after the 100th step of a QW with |Aq;| = 2, the 
pulse duration must be shorter than T^g gg(200) = 0.21 ns. 

However, applying the photon kick when the ion is at its turning point, i.e. the ion 
is slowest and a is real, the scaling is less demanding. We find 

Tf=o.9g{\a\) = exp (-17.03 - 0.021n(|a|) - 0.1 (In (|a|))^) . (36) 

Most importantly, the prefactor of the term linear in ln(|a|) is much smaller than for 
an imaginary a. For the 100th step, the pulse duration therefore only has to be shorter 

§ Equations ([35|l and (|36p are the results of quadratic fits to double- logarithmic plots of pairs (Tp, |q;|) 
for a fidelity / = 0.99 and |a| < 10. Higher motional amplitudes were not considered due to sizable 
additional numerical effort. For the following estimates, the scaling is considered to be preserved. 
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than TJLq gg(200) = 2.18 ns, which is within the specifications of a fast-switching electro- 
optic modulator and our current continuous-wave laser system. Timing the application 
of the photon kick to the (spatial) turning points of all the coherent oscillations occuring 
during the QW is possible, because we start the QW in the motional ground state and 
the position states are aligned along a line in the co-rotating phase space. Thus, coherent 
states of different |Aa| reach their turning points simultaneously. 

Given our experimental parameters, in particular the width of the ground state wave 
function zq = 10 nm, the coherent motional state of maximal amplitude, \amax = 200), 
would have a real-space amplitude of {amax \ z la^ax) = 4 ]im. At such high motional 
amplitudes anharmonicities of the trapping potential must be considered. These depend 
on the design of the electrodes and could be eliminated, e.g. by designing the Paul trap 
electrodes in a hyperbolic shape [5lj. Additionally, micromotion might increase the 
deviation from the ideal walk, for example by reducing the overlap of the additionally 
oscillating motional wave functions. However, it will remain negligible when the QW is 
implemented in the axial degree of freedom of an ion in a linear Paul trap. 

6.2. QW in higher dimensions 

A QW in two or three dimensions is possible by additionally considering the motion 
in the radial direction. The pulse sequence for a step of a QW is then the subsequent 
application of the shift operator in each direction where each operation is preceeded by 
a coin toss. 

More possibilities and reduced technical requirements might be achieved by 
trapping more than one ion and considering the collective degrees of motion in one 
direction. Reference [52] describes the scheme with two ions, creating a 4-sided coin, 
where two coin states affect the walk in the center-of-mass motional mode and the 
other two in the stretch mode of motion. In particular, possibilities with the coin being 
initialized in an entangled state are investigated. 

A photon kick, as decribed above, induces motion in all motional modes in the 
direction of the effective wave vector k, according to the respective coin states. That 
is, with ions, one step of the QW consists of a coin operation on the 2^-sided coin 
and a single shift operation, that displaces the part of the motional wave function 
related to each coin state into opposite direction in phase space of one motional mode. 
The particular difficulty is to assign the 2^~^ pairs of coin states to A^ different axial 
motional modes such that for each coin state the corresponding state dependent force 
induces motion in one direction in one certain mode [53]. One possibility to obtain 
the required number of motional modes is to add ions that do not contain a transition 
corresponding to the coin states and therefore are not affected by the photon kicks. 
That would be in our case ^^Mg ions without hyperfine structure and therefore no coin 
states and no corrsponding transition. With this, the implementation of a QW in four 
dimensions is possible using three ^^Mg ions and one ^^Mg ion. For more dimensions 
the issue arises that the ^"^Mg and ^^Mg ions have to be arranged in such a way that for 
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each coin state the corresponding state dependent force induces motion in one certain 
mode, which has not been clarified yet. 
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